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1. Introduction 

Is "staggered QCD" really QCD, or is it just a model of QCD ? This is the question addressed 
in this note, and it is worth emphasizing two points right now: the issue is whether the staggered 
approach yields the right continuum limit and the question is asked at the non-perturbative level. 

The staggered action is not doubler-free; the Kogut-Susskind construction manages to reduce 
the fermion multiplicity from 16 to 4 [[Tj]. To simulate QCD with Nf - = 2+1 dynamical quarks 
staggered studies employ the "rooting trick", i.e. the Boltzmann weight used in the simulation is 

e S af = det l/2( Dstffl J A ^l\ DAms ) e So (U) 

where m uc j = (m u +m,i) /2 is an average up and down quark mass, m s is kept at the physical strange 
quark mass and Sq denotes the gauge action. Throughout, I use the abbreviation D st m =D st +m for 
the massive staggered Dirac operator. The "rooting trick" is the way how the unphysical "tastes" 
(see below) are excised from the sea sector of the theory, and the issue of debate is whether the 
effective action may be re-expressed in terms of an artificial undoubled formulation which is local. 
Ideally, this construction would hold exactly on typical backgrounds, but in view of the first point 
emphasized above it might be OK if one can prove that a slight mismatch in the determinants is 
just a cut-off effect. I will explain below why it is so useful to have a local underlying theory. On 
the level of the effective action the question is whether multi-plaquette correlators have unphysical 
cuts in the limit where the lattice spacing a is sent to zero. 

The staggered action not being doubler-free has an impact on the valence sector, too. Here, 
sophisticated reduction techniques have been developed (already in the quenched era) to project 
out the desired spinorgtaste combination from a given staggered operator. The problem is that the 
way how the unphysical tastes are excised from the valence sector is (in general) slightly different 
from the rooting procedure employed in the sea sector. Evidently, this means that at finite lattice 
spacing the theory is not unitary. Again, this does not necessarily lead to a disaster, if one can show 
that the non-unitarity is a genuine cut-off effect. 

Personally, I tend to see the locality as the crucial issue and the unitarity as a point of lesser 
severity. This is so, because there are reasons to believe that a solution to the first problem would 
automatically provide a bypass to the second one, as sketched below. 

The controversy, as I see it, is thus a debate whether the rooting procedure in the sea sector 
would introduce spurious degrees of freedom which corrupt, at the non-perturbative level, the long- 
range properties of the effective action in such a way that something survives the continuum limit. 
As of now, the answer is not known. All I can provide is an attempt to explain basic facts needed 
to understand the origin of the problem and then collect, in a second step, pieces of evidence in 
favor and against the legality of the staggered approach. It turns out that the evidence "in favor" 
outnumbers the one "against" by far, but in the end the sheer number has no say. 

A key issue in practical applications of staggered fermions is to offset the taste symmetry 
breaking as much as possible. In the course of the simulation this is done by using improved/filtered 
varieties of the staggered action which may reduce the taste splitting quite drastically (see below). 
In the analysis of pseudo-scalar meson properties Staggered Chiral Perturbation Theory (SXPT) is 
used, an effective field theory that allows one to parameterize and thus "take out" the leading taste 
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breaking effects. Many of the pieces of evidence to be presented below will use one or the other 
strategy to tame the taste symmetry breaking. Still, from a fundamental viewpoint it is clear that 
either one of these techniques is immaterial - if the staggered approach is correct, it will yield the 
right continuum limit for arbitrary observables without any of these. 

Obviously, there are two ways how the controversy could play out. One possible outcome 
is realized, if someone manages to give an unambiguous (mathematical) proof that QCD with 
Nf = 2+1 staggered quarks is in the right universality class. The other option is, of course, that 
someone manages to find a single observable where he/she can demonstrate that the staggered 
answer, after a continuum extrapolation, is wrong. 

Before getting into details, let me comment on one trap to be avoided. It is true that some of the 
more recent studies with Nf = 2+1 dynamical staggered quarks agree (in particular in observables 
relating to T spectroscopy) with experiment to a precision that is currently not matched by another 
fermion discretization (see e.g. [§]). Still, from a conceptual viewpoint this is no help in deciding 
whether the staggered approach is fundamentally correct or "epicycle physics". In particular the 
agreement of a mixed QCD/electroweak observable with experiment is not really an argument, 
since what we want to find out is whether there is a contribution from beyond-the-SM physics in 
such observables. What comparing to experiment is supposed to tell us is whether the Lagrangian 
relevant to a problem is the one we stick into our code, and for this the method itself has to be from 
first principles. Hence, whether QCD with Nf = 2+1 staggered fields is fundamentally correct is 
an issue that needs to be solved theoretically - by pure thought or with the help of CPU power. 

The plan is to first review some elements of the staggered action and taste representation. 
The question about the legality of the "rooting trick" is then exposed as a locality problem that 
would be solved if one could find a "candidate" operator that satisfies certain criteria. In the free 
theory the search has ended successfully, and I discuss the four constructions known to date. For 
the interacting theory I assess pieces of evidence that come from studying spectral and eigenmode 
properties of the staggered action in 4D. For many observables in 2D it is easy to reach good 
statistical precision and this leads to a somewhat more mixed picture in the Schwinger model. Still, 
just comparing the staggered and the overlap contribution to the effective action [i.e. correlating 
— ^logdet(D st ; „) and — logdetD ov m in a scatter plot for ra>0] lets one ask whether D ov m might 
actually represent (one version of) the desired "candidate" action. Finally, a few words on what 
can be learned from analyzing data with Staggered Chiral Perturbation Theory are in place. The 
problem being unsolved, I refrain from drawing any conclusion and end with a summary. 

2. Review: Staggered action and taste representation 

The staggered action derives from the naive fermion discretization 

a 3 

S na = —Y,Y(x)Yti[U t i(x)Y(x+fi)-Ul l {x-tl)Y^-^)] (2-1) 

which is known to lead to 4 (in 2D) or 16 (in 4D) fermions in the continuum (2 d species in d 
dimensions). The observation by Kogut and Susskind was that the x-dependent field transformation 
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\lf(x) = y(x)x(x),Y( x ) = X{ x )y{ x )\y{ x ) = i^iii^it withjc=raa and the staggered phase factor 
Tj Ai (x) = (-1)£m<v"v brings it to the form 



a 3 



^na/st = y I^WinW [^W^( x +A)-^(^-A)z(^-A)] (2-2) 



where the 2 (in 2D) or 4 (in 4D) components of % (in general: 2 rf / 2 ) manifestly decouple. Therefore, 
by just "downgrading" # to be a 1-component field, the action (^2|) is changed to the staggered 
one which yields only N t = 2 d l 2 fermions in the continuum. The Kogut Susskind procedure thus 
manages to "thin out" degrees of freedom, albeit at the price of distributing/intertwining "spinor" 
and "taste" (the latter being the modern word to address the full fermion content of a staggered 
field). An important point is that the SU (N t )A taste symmetry in the formal continuum limit reduces 
to the abelian invariance %(x) — > e ldAr,5 ^x(x),%(x) — > e~ l9Ar,5 ^x( x ) of the massless action at finite 
lattice spacing, and this is good enough to prevent additive mass renormalization. However, further 
exact "thinning" by purely algebraic means is impossible, since the Kogut Susskind procedure has 
exhausted the eigenvalue degeneracy of the naive action (the latter has only 2 rf / 2 -fold degenerate 
eigenvalues on typical backgrounds, not 2 rf -fold, as one might have naively guessed). 

Modern staggered simulations use an individual staggered field for each continuum flavor to be 
studied, i.e. QCD with Nf = 2+1 quarks requires the ^-fields u,d,s. The taste structure of physical 
operators is non-trivial ^ |J]. The underlying taste identification may be attempted in momentum 
space pp or in coordinate space H 0, §1. The latter version builds on the hypercubic decomposi- 

, A A A A, J 

tion x(x,x+a\, ...,x+a\+a2+a3+a4) — > q(X) which collects all spinor/taste information in a 2 
sublattice (of the original one) into a "masternode", i.e. a point in the blocked lattice with spacing 
b = 2a. Denoting these coarse points with uppercase symbols, the free action reads 

Sst = b 4 £ q(X) [V M (7 M ®/) - ^A^Ysm^s)} q{X) (2.3) 
x,n 1 

where q{X) is a 16-component "quark field" (in 4D) and the first and second derivative operate on 
the blocked lattice, i.e. (V^)(X) = ^ x +m-<i(x-m and ( A ^)(X) = gg±Mlzjffl±gg^ _ The 

tensor structure in ( |2.3[ ) refers to spinor/taste, where Tu = y^fo = J$, i.e. another set of gamma 
matrices is used to map out the taste algebra. The main observation is that in this basis the trouble- 
some taste changing interactions stem from a dimension 5 Wilson-type term which is "skewed" in 
spinor/taste space (and hence not positive semidefinite). 

The case of interest is the interacting theory, and here the question is whether all taste changing 
effects would go away, in the continuum limit, without leaving any trace. Of course, from a formal 
viewpoint this seems almost obvious, since it is a dimension 5 piece 1 in ( ^3| ) that mediates the 



taste changes, and irrelevant operators have no anomalous dimensions. As pointed out in [jlOp the 
whole issue is a peculiar exchange-of-limits question. If the formal continuum limit is justified 
(and the view that naturally derives from the taste basis thus adequate), then everything is fine, 
indeed. The problem is just that one would need to prove first that the staggered action yields the 
correct continuum limit to be sure that adopting this viewpoint is not misleading. 



As worked out by Luo the 0(a) term in (2.3) can be shifted to 0(a ) by a change of basis | 
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Figure 1: Staggered and overlap eigenvalues on 4 selected configurations in 2D. On benevolent config- 
urations the improved taste symmetry achieved through filtering results in a 2-fold (4-fold in 4D) near- 
degeneracy and the right number of near-zero modes. To aid comparison the low-energy overlap spectrum 
has been mapped onto the imaginary axis with a stereographic projection. Figure from [|39||. 



One way to think about taste changing interactions is to notice that in the interacting theory in- 
dividual taste (or spinor) components of a g-field still sit out there in the little 2 4 hypercube attached 
to a blocked "masternode". They see slightly different local gauge fields (therefore spinor(g)taste 
projection operators must be gauged, e.g. the staggered 75 is a 4-link operator). Equivalently, in the 
momentum space view each spinor/taste component has a reduced Brillouin zone with extension 
2%jb = 71 /a in each direction and this means that hard gluons with one component ~ ft/a will flip 
one taste to another one. In any case it is clear that the identification of physical and staggered 
flavor may only work, if taste interactions are minimized (ideally: eliminated). 

Everything below will hinge on this taste symmetry (which at finite lattice spacing is only 
approximate), and Fig. |l| is an attempt to convey a down-to-earth view. If the taste symmetry were 
exact the staggered eigenvalue spectrum would show a 2-fold (in 2D) or 4-fold (in 4D) degeneracy. 
On typical configurations smearing (i.e. using a filtered fermion operator, see below) is a convenient 
tool (read: cheaper than decreasing a) to achieve a rather good A^-fold near-degeneracy, and this is 
the basis of the rooting procedure. In other words, if the symmetry/degeneracy were exact, it would 
be conceivable that one might reduce the degrees of freedom by another factor of 2 d l 2 (analogous 
to Kogut and Susskind, but with different technical means). However, at finite lattice spacing there 
are always some configurations where the taste symmetry is not approximately restored (rightmost 
panel), and the issue concerning the "rooting trick" is basically the question whether these "bad 
guys" would retain non-zero measure in the continuum. 
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3. Problem: rooting versus locality 

The earliest warning regarding the legality of the "rooting trick" that I am aware of came 
in a paper by Marinari, Parisi and Rebbi in 1981: "On the lattice the action 5g — |trlog(D st m ) 
will produce a violation of fundamental axioms, but we expect the violation to disappear in the 
continuum limit and then recover the theory with a single fermion" tjiTp. 

As mentioned in the introduction, the main problem is the potential loss of locality at finite 
lattice spacing. With any undoubled Dirac operator, e.g. D=Dw, one has 

Z\ Nf = Jd[U,yM e- s °W-^* D v = Jd[U] det N f{D) (3.1) 

where the first sum extends over Nf flavors, i.e. an integer power of the determinant emerges rather 
naturally with Nf degenerate fermions. The question is whether one may take fractional powers if 
a formulation with doublers yields too many fermions in the continuum, specifically whether 

Z\ Nf = Jd[U] det^/ 4 (D st , m ) e- s °M (3.2) 

is a legal implementation if Nf is not a multiple of 4. The point is that it is not clear whether the 
resulting effective action S e s = Sq— -j-trlog(D st)/n ) receives contributions from unphysical degrees 
of freedom which persist in the continuum. In principle, such a disaster might happen with any 
non-zero prefactor, and locality is thus a issue in the staggered approach for any Nf > 0. 

The reason why locality is so important is universality. We are all used to the notion that we 
do not need to worry about the details of the action used, since the continuum limit is a 2 nd order 
phase transition in the course of which all physical correlation lengths diverge and any knowledge 
about irrelevant pieces in the action is washed out. However, the standard theorems ensuring these 



nice properties build on the fundamental action being local (at finite lattice spacing) Q12 1 . 

The soundness of the effective action would be guaranteed if one could "integrate in" a fermion 
in (3.2), i.e. if one could find a legal doubler-free "candidate" action D ca m with [ |l3| ] 



det 1 / 4 ^^) = JDWM e-/<^W . ( 3 .3) 
The attribute "legal" means that it satisfies the requirements 

Dc, = i7nP^ + 0(p 2 ) (3.4) 

I |A»(*,y) 1 1 < C<T v|v ~ v|/o with C,v independent of U (3.5) 

\\SD ca (x, yi U)/SU(z)\\ < Ce- vmin (M'l^l>/ a (ditto) (3.6) 



where (3.5 3.6) are designed to ensure locality in the space of couplings and in the gauge field, and 
"doubler-free" is formalized by demanding that D ca {p) is invertible at all non-zero momenta (mod 
111 /a). Of course, one might allow for cut-off effects in (fO|), i.e. one might relax it to 

det(D cai , n ) ^const-det^Dstm) (3.7) 
but the problem will be to prove that the difference is really just a cut-off effect. 
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Figure 2: Localization function f(r) in physical units versus r/ro for the "candidate" action of Jl5|]. The 
local logarithmic derivative (the effective Hoc) scales and the operator is thus non-local. Figure from @. 



What might be a promising "candidate" action ? Sometimes it is easier and more instructive to 
discuss first what is not a useful strategy. The staggered "rooting problem" is of this type, since it 
is evident that simply taking the fourth root of D stm is not a viable option. Such an operator would 

' 1/4 

have the right determinant, but it would violate other requirements, in particular D ca = D st with 
D st the massless staggered operator would fail to obey (3.4). 

That the technique [14] is available to decide whether a potential/proposed "candidate" action 
obeys the locality bound ( fO| ) has bee shown in [JTJJ]. These authors consider a slightly different 
"candidate" operator, D cam = {Dl tm D st m )H 2 = (—D 2 t +m 2 )U 2 > 0, and measure the localization 
function f(r) = sup{ HvMlb | y\\i = r } where y(x) = £D(jc,z)il(z) with Tj a normalized 
random vector aty. For a local D one has f(r)<xe~ r ' noa , ideally with r\ oc °<a, but any a-dependence 
with ^phys/loc — >oa is fine. Employing matched lattices they find that /(r/ro) scales (and there is no 
convincing exponential fall-off pattern, see Fig. ||), i.e. the operator is non-local. Thus, apart from 
not qualifying as a fermion action since its spectrum lies on the real axis, this "candidate" operator 
may also be ruled out for more sophisticated reasons. Using an improved/filtered staggered action 
in the above-mentioned construct does not alter the conclusion [16], which is hardly surprising, 
since in the previous work the continuum limit has already been taken. 

Still, the question remains: is there any "candidate" operator satisfying all requirements ? Note 
that finding a legal 1 -flavor D ca „, is sufficient to show that staggered fermions yield an acceptable 
effective action. The continuum limit might be OK even if this is not the case, but in such an event 
a proof that the approach is safe would - almost for sure - require truly novel concepts. 



4. Free case: four constructions 



Given the importance of the issue whether a "candidate" action exists that satisfies all relevant 
criteria, it is encouraging to see that at least in the free theory the answer is positive and the task to 
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construct such an operator has been completed. Interestingly, the solution is non-unique. 

In [|T^] Adams starts from the observation that the massive staggered operator in the taste basis 



( |2.3| ) on the blocked lattice may be used to build the combination 



D' D 



V 2 + ^-A 2 + m 2 



(4.1) 



which is simultaneously diagonal in spinorigitaste. Analogously, a free generalized Wilson operator 



D, 



^n7n + |W + m on the blocked lattice yields 



D J D 



-W + m) 



without taste structure. With det(Dj tm ) = det(D stm ) and det(D f 



■ det(D 



W,m/ 



(4.2) 

it follows that 



D C a,m would trivially have the right determinant, if the two square brackets in (|4JJ) and ( |4.2| ) agree 



— A 2 + m 2 = ( T U' +/;/ 



det 1 / 4 (D stjm ) = det( J D cai „ 



(4.3) 



Hence, a natural "candidate" operator in the free theory is D ca = V^y^ + y ^ A 2 + m 2 which lives 
on the blocked lattice. The less trivial part is to verify its locality. Using a technique developed by 
Fujikawa and Ishibashi he shows that the couplings fall off exponentially with n oc = y/Sa/m 
and this means that for m > the operator is local. 



a 10 







In [ |18P Maresca and Peardon investigate the massless case. They start from a "candidate" 
operator D ca = iy^P^ + Q on the blocked lattice with local (i.e. not ultra-local) P^ and Q. Again, 
det(D ca ) = det 1 /^' (D st ) is trivially realized if they can ensure that Dl a D c 



aim they consider the ansatz P^ = Lr>oL|rf|<r Mpjifay) an d ditto for Q in terms of a tower of 
operators with ever larger footprint, i.e. CBp?n,COq' have range r. Ensuring the main criterion does 
not uniquely pin down all the coefficients, and they use this freedom to optimize D ca for locality. 
This optimization is performed numerically, both in 2D (A^ = 2) and in 4D (A^ = 4). With this 
strategy they make sure that D ca has the right determinant [my understanding: up to the exact 
(non-topological) zero-modes], and it turns out that the resulting operator is indeed exponentially 
localized. Presumably triggered by the spectrum of their "candidate" operator in 2D (see Fig. 0) 



» rtaste 

>> l N t 



D' st D, t . To that 



they went on to ask whether they could make it obey the Ginsparg Wilson relation [19]. Indeed, 
considering only Q of the form Q = —/SR with a local R makes the "candidate" operator satisfy 



{D ca ,Y 5 }=D ca 2Ry 5 D c 
but the price to pay is that the localization range increases. 



(4.4) 



In Q20| ] Shamir uses the renormalization group (RG) framework to deal with the free case. The 
general idea is that even in the interacting theory performing n blocking steps (for simplicity: n 
coarse grainings by a factor 2) creates a "two cut-off" situation where the physics accessible on 
the coarse lattice is given by the final cut-off cT x /2 n ' , but the taste splitting is still given by the 
original aT x . Thus, starting from an ever finer lattice and applying more blocking steps is a way 
to improve the taste symmetry (cf. talk by F. Maresca [|2~T|]). Infinitely many steps render it exact, 
i.e. D n -^Doo®h, while D ca = D n after n steps satisfies det 1 / 4 (D st ) = det(D ca )det 1 / 4 (r) where T 
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Figure 3: Eigenvalue spectra of two "candidate" operators D ca on the blocked lattice in 2D which reproduce 
det 1//2 (D st ) of the massless staggered action in the free theory. On the left D ca is exclusively optimized for 
locality, on the right the (standard) Ginsparg Wilson relation is imposed as a constraint. Figure from | 



contains only (original) cut-off excitations and should maintain the Symanzik class, i.e. det(T) = 
const-(l+0(a 2 )). The massless staggered action on the original lattice satisfies {Do, (Ys'SiTs)} = 0. 
Interestingly, the action after n blocking steps satisfies a generalized Ginsparg Wilson relation [ |20| ] 



{D n ,(Y 5 ®T 5 )} = D n —(y 5 (g)i; 5 )D n 



(4.5) 



where a n is a parameter of the blocking procedure. If ( fO| ) were complemented with a (75(8175)- 
hermiticity of D n , one would easily obtain D n +D] 1 =D n ^-D] v i.e. that the spectrum of D n lies on 
a Ginsparg-Wilson circle. However, it is straightforward to derive spectral properties without such 
a (y5®T5)-hermiticity. Start from the representation D n = £y Ay^-Vy in terms of bi-orthogonal R/L- 
eigenvectors, i.e. v\u k = 5,-^- Upon sandwiching the GW relation (4^5) with vj and u k , one obtains 
XiV i (...)ui ( + vJ(...)ui ( Xk = Xiv\^r(...)iikhk where ... stands for 75(8)15. Thus, for an arbitrary pair 
(i,k) of eigenmodes, one has vj(7 5 ®x 5 )u k =0 or Xi+Xu— ^-A,-A^=0. In particular for i=k it follows 
that Vy(75<8T5)w ; - 7^0 necessarily implies 2Xj — -^-Xj = 0. Hence the spectrum crosses the real axis 
at Xj € {0, a n }. The more general situation with a spectrum between two "Hasenfratz-Niedermayer 
circles" (tangent to the imaginary axis) is avoided, and the reason why is simply that Shamir uses 
a non-overlapping blocking kernel. 



There are two more articles to be mentioned. In Q22Q Giedt gives an exploratory discussion of 
the interacting case, based on RG concepts. In [^3|] Neuberger uses a different setup; he recasts the 
question of the validity of the "rooting trick" into a local field theoretical framework in 6D. 



5. Interacting theory: spec(D st ) in 4D 

Before discussing what can be learned from analyzing spectral properties of the staggered 
Dirac operator, let me recall the concept of filtering. What is meant is simply a change of the 



covariant derivative in the fermion action [ ]24| , |25J, pq, 27 



8|, [29p through a "fat link" recipe like 
UYP 'x) ¥ (x+fi)- ¥ (x). (5.1) 
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Figure 4: Chirality versus eigenvalue of staggered eigenmodes. On coarse lattices there is no distinctive 
signature, but on somewhat finer lattices rather aggressive filtering reveals continuum like features. Near- 



zero modes are predominantly chiral, non-zero modes are not. Figure from [33] 



The idea is that such a modification alters/diminishes cut-off effects without changing the Symanzik 
class, i.e. a theory with 0(a 2 ) cut-off effects will retain this characteristics, but the extrapolation 
slope is changed (typically reduced). For staggered quarks the obvious design goal of a filtering 
prescription is to render D st "immune" against p^n/a gluons, and the question is how efficient a 
given recipe is at ameliorating the taste symmetry. For Wilson fermions the additive mass renor- 
malization gets reduced through the modification ( 5A ) and for related technical reasons even the 
overlap operator [30] benefits from such a filtering [ ]3"T| , |32| ]. 

An excellent starting point is found in a talk by A. Hasenfratz [33]. In the (formal) continuum 
D\j/=X\l/ Dj5Y = ~ ^JsY, tnus non-zero modes come in complex conjugate pairs (A is purely 
imaginary) and have zero chirality (^ t y 5 i/A = if X ^0), while zero modes have non-vanishing 
chirality (£^£ = ±1 if D£ = 0). On coarse lattices unfiltered and modestly filtered staggered 
quarks show no indication of such "continuum like" behavior (see left panel of Fig. |J), but with 
cT l ~2GeV and some rather aggressive filtering the situation changes, on a qualitative level at 
least (see right panel of Fig. ||). The idea is that the same picture would emerge with the "thin link" 
action, too, if one could push to much smaller lattice spacings (say cT l ~20GeV or so). At least it 
becomes evident that cut-off effects may completely mask the underlying continuum theory. 

A more complete investigation has been carried out and presented last year by Follana et 
al. [34, [$5|]. The difference is that they show results for individual configurations and compare 
several filtering recipes and gauge actions (at fixed lattice spacing). Furthermore, they determine 
the (gluonic) topological charge q of their configurations. With lattice spacings a~ 1 ~ 2 GeV only 
the combination of an improved gauge action and a highly filtered D st allows to see a separation 
into near-zero modes with non- vanishing chirality and non-zero modes without chirality (cf. Fig.|5|). 
Furthermore, the number of near-zero modes is 4\q\ (i.e. 2\q\ on either side of the spectrum), and 
this means that staggered quarks satisfy an approximate index theorem. 

The next step has been taken in a paper together with Christian Hoelbling and Urs Wenger 
where two new features come in [31]. The first one is that we compare to the overlap spectrum on 
the same configuration. Apart from finding 4\q\ staggered near-zero modes on typical backgrounds 
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Figure 5: Chirality and eigenvalue versus "serial number" of the lowest few staggered eigenmodes of the 
positive staggered half-spectrum. The left panel is for Wilson glue, the right one for the Symanzik action, 
with a~0.1 fm in either case. The staggered action with "HISQ" filtering sees \q\ =2. Figure from [53], 
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Figure 6: Same as the lower part of Fig. 5, but with two different lattice resolutions (keeping the physical 
volume fixed) and with overlap eigenvalues on the same configuration for comparison. After rescaling with 



Z^'/Z^ v one sees a quantitative single-overlap to staggered-quartet agreement. Figure from [bin. 



with overlap charge q we show that after rescaling with the relative Z^-factor one finds a quantitative 
agreement of staggered quartets with the (non-degenerate) overlap eigenvalues (cf. Fig^). The 
second new ingredient is that we compare several matched lattices (with fixed physical volume) 
and investigate how the intra-quartet splitting decreases with a finer resolution (see below). 
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Figure 7: Cumulative eigenvalue distribution for two filtered D st and the "thin-link" D m . Figure from [pq]. 
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Figure 9: Same as Fig. 8, but separated according to \q\ and with much better statistics. Figure from [37 1 



For quite some time it was believed that staggered quarks are insensitive to topology and would 
never reproduce the sectoral random matrix theory (RMT) predictions for the distribution/ratio of 
the lowest eigenvalues. Given the correspondence to the overlap spectrum it is clear that this cannot 
be true. In [36] it is shown that on sufficiently fine lattices and with filtering the agreement is just as 
good as with the overlap action (see Fig. [7]). In [Q] it was already shown that the ratio of low-lying 
eigenvalues [i.e. (A,-) / (A/) with A, denoting a staggered quartet mean] satisfies the RMT prediction, 
up to small violations which might be finite volume effects or remnant discretization errors. 

The last point to be covered is the evidence that the staggered taste symmetry violation might 
indeed be an 0{a 2 ) cut-off effect. A first attempt was made in [|1|] where the taste splitting is de- 
fined through the "spread" of a staggered quartet that is to mimic a (single) overlap mode. Plotting 
the splitting versus a 2 , everything seems consistent with the hypothesis that it might disappear in 
proportion to a 2 [with filtering, without one is certainly not in the scaling regime], see Fig.|8[ In 
[37] these findings have been corroborated with much better statistics (see Fig.|9|) and in fl2"ij ] a 
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more physical definition of the taste splitting has been found to lead to the same overall picture. 

In summary, investigations in 4D have so far revealed evidence in favor of the "rooting trick". 
The problem is just that even the most convincing numerical evidence does not amount to a proof. 

6. Interacting theory: ^ sca ,Zto P in 2D 

The main reason for working in 2D is that for a large number of observables a rather impressive 
statistical precision can be reached. Hence, if there is a problem with the staggered approach, one 
might think that it is most likely to be detected in such a setup. 

In 2D a square root of the determinant is used for an odd number of dynamical fermions, since 
the staggered action yields 2 fermions in the continuum. Results below are in the (generalized) 



Schwinger model [38], i.e. QED in 2D with Nf massive fermions. Note that the fundamental 
charge g has the dimension of a mass and does not run. Therefore, g may be used to set the scale, 
and in this case the relation j8 = 1 / (ag) 2 holds exact. For Nf = 1 there is the analytic result [ [38| ] 

lim^-^ = — = 0.1599... (6.1) 

mjo g In*! 1 

which is due to the axial anomaly and does not signal spontaneous symmetry breaking. 

The challenge is to see whether one may directly reproduce ( |6.1[ ) with staggered fermions. To 
this aim we sample quenched, compute all eigenvalues (both of D st and D ov to have a cross check, 
and in either case with several filterings) and determine the observable in the analysis program via 



(6.2) 



g V (det^D-)) ' g 2V (det"// 2 (D£)) 

vyt 

det(D°:)=l\((l--)X+m) , det(D^)=n(A+m) (6.3) 

where A = (1/A — 1/2) denotes the stereographically mapped overlap eigenvalues. The pref actor 
A in the staggered expression is designed to compensate for the 2 tastes contributing; this is the 



valence counterpart to the "rooting trick" in the sea sector. The bare condensates (6.2) have a 



logarithmic divergence of type m/g-\og{ma), i.e. the free case needs to be subtracted, but there 



is no multiplicative renormalization. For technical details see [39, 40]. All plots below will have 
»?sea = 'Wvab i.e. at least with D=D 0V the theory is exactly unitary (at finite lattice spacing). 

Given this setup it is clear that we can draw "continuous" la (statistical) error bands for the 
subtracted condensate xLa/s as a function of m/g. The volume is chosen such that the Compton 
wavelength of the lightest degree of freedom in the chiral limit for Nf = 1 fits 4 — 5 times into the 
box. The first observation in Fig. [l0| is that our overlap data satisfy the continuum expectation 

Xsca/ S x * 



g/ m (N f =0) 

const (N f = l) (6.4) 

m/g (N f = 2) 



at least on a qualitative level. Moreover, the overlap results with and without filtering look almost 



identical. And finally, in the chiral limit the 1 -flavor result is consistent with the exact value (6.1). 



With staggered fermions the picture looks rather different. The most obvious problem is that in 
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Figure 10: Scalar condensate versus quark mass (everything in physical units) in the Schwinger model, with 
D s t and D ov , with Nf = 0, 1 , 2, without (left) or with filtering (right). Figure (with modification) from 



the chiral limit these condensates tend to zero for all Nf, in marked contrast to (6.4). However, 



at the coupling used a single filtering step suffices to dramatically improve the situation - for not- 
too-small quark masses at least. Still one has lim m _>o/£s C a/£ = at finite lattice spacing (even for 
Nf = where the slope is so steep that already the first pixel is off scale), but the point where this 
happens has been shifted to the left. This is in line with our previous observation that the splitting 
within staggered eigenvalue doublets gets dramatically reduced through filtering (cf. Fig. [I] and 



Ref. [p9|]). However, there is a large numerical factor relating the intra-taste eigenvalue splitting and 
the breakdown mass, and this ought to be understood before one claims a relationship. The main 
question that emerges from the r.h.s. of Fig. [R] is whether one can reproduce ( |6.1[ ) with staggered 
fermions if one first takes the continuum limit (at finite m/g), pushing to the chiral limit in a second 
step. It turns out that this is indeed the case, viz. 

y mlim *SWlS) = JL , UmUm jgWg^> =0 ( 6 .5) 

H.H^rW^l.jL , limlim O"fe° 2 ) = 4_ (6.6) 

m^Oa^O g 2K i > 2 m^Ofl^O g 111 5 ' 1 

where every identity except the one with the exact zero is to be read as "our data are consistent 



with the hypothesis that ... " [40]. Thus there is a non-commutativity phenomenon specific to the 
staggered action. In other words, for this theory the staggered action is not in the right universality 
class, if one implements Xscn/g through ( |6.2| ) and insists on a particular "bad" order of limits. 
Perhaps more surprising to staggered skeptics is that the "good" order of limits (first a — ► 0, then 
m — > 0) brings it to the right universality class. This means that staggered quarks - with the help of 



an extra guiding principle - do perceive the chiral anomaly Q40|]. 

What makes staggered quarks fail, in this specific case, with the "bad" order of limits ? We 
did two experiments. Already in [ |39| ] we asked what happens if one combines staggered sea with 
overlap valence quarks or vice versa (BTW, [|^] contains the first published "hybrid action" test). 
In such a setup one has the "taste reduction" problem only in one of the two sectors. We restrict 



ourselves to m sea =m va i. As shown in Fig. 11 the results with overlap valence quarks on staggered 



sea look similar to the genuine staggered data shown before. Again, the curves tend to zero below 
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Figure 11: Same as Fig. 10, but for the formulations overlap valence on staggered sea (dark) and staggered 
valence on overlap sea (light), without (left) and with filtering (right). Figure (with modification) from 
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Figure 12: Attempt to "zoom in" on the coefficient of the divergence in the — Q case. Overlap fermions 
yield the right result, regardless whether m — > or a — > is taken first. With staggered fermions another 
non-commutativity shows up which does not involve a rooted determinant. Figure previously unpublished. 



some "critical" quark mass, and the latter happens to be more pronounced (and smaller) with filter- 
ing than without. An interesting observation concerns the reverse combination, staggered valence 
quarks on overlap sea (which in 4D would be a crazy thing to do). Here, the curves blow up at 
the same "critical" mass. At least in this example the notion of exact chiral symmetry in either 
the valence or the sea sector alone does not enable one to go closer to the chiral limit (at fixed a) 
than with the generic "staggered on staggered" combination. This is potentially bad news for mod- 
ern attempts to evaluate overlap/domain-wall valence quarks on "cheap" staggered ensembles (see 
pi) , |42| ] for a discussion of specific issues and [43, |4|] for the effective theory). For our problem 



the bottom line is that the failure of the staggered action in the 1 -flavor theory with the "bad" order 
of limits cannot be attributed to the sea or valence sector alone. 

The second experiment is the result of a discussion with Claude Bernard. To find out whether 
the non-commutativity (5.5,6.6) is causally connected to the "rooting trick" in the 1 -flavor theory 
with staggered fermions, we checked whether one can compute the coefficient of the quenched 
divergence in (6.4) with staggered fermions. In Fig. 12 the product mXs^° V is plotted versus m/g. 
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Figure 13: Topological susceptibility Xtop/g 2 m the 1-flavor theory with det 1 / 2 (D stim ) (left) or det(D ov m ) 
(right) at fixed m/g= 0.1, plotted versus (ag) 2 . The continuum limit seems to be universal. Figure from 



With overlap quarks the order of the chiral and the continuum limit is irrelevant. With staggered 
quarks taking wi— >0 first yields a vanishing result, hence taking a—>0 in the second step leads to an 
exact zero. However, the result is different if one reverses the order of limits. With a— >0 the "dip" 



in Fig. |12| gets filled in for every finite m/g. Therefore, taking m— >0 in the second step one ends 
up with a finite result (which is conjectured to be identical to the coefficient obtained from overlap 
fermions in the combined chiral and continuum limit). 

At the moment it is not clear what the precise implication of any of these experiments is for the 
case of interest (QCD with Nf =2+1 quarks at their physical masses). Evidently, the chiral limit at 
finite lattice spacing seems to be troublesome for staggered fermions. This has already been seen 



by Smit and Vink in Ref. [45] where they determine (mass dependent) renormalization factors for 
the topological charge as seen by the massive D st m . These factors tend to diverge for m— >0. Our 
contribution was just to promote this into troubles with a real physical observable. In any case the 



bottom line is that the non-commutativity (^5^6) is not tied to the "rooting trick" per se but rather 



reflects a general mismatch between the sea and valence sector. This fits the view by Bernard who 
observes that similar non-commutativities show up in the effective theory for staggered fermions 
in 4D with partial quenching, i.e. if the sea and valence quarks have different masses [fig]. 

Finally, I should add that the non-commutativity ( |6.5| , |6T6| ) is more the exception than the rule. 
We checked a variety of other observables in the Schwinger model and did not find anything similar. 
To give an idea Fig. [13| shows our data for the Nf = 1 topological susceptibility Xiop/g 2 at fixed 
physical quark mass. Having data with three filtering levels and with several couplings in the 
Symanzik scaling regime allows for a fairly safe extrapolation to the continuum. The results with 
the two discretizations seem consistent, and also the exact continuum result lim m ^oZtop/g =0 is 
reproduced with staggered fermions, regardless whether m — > or a — > is taken first ||40|]. 

7. Correlation of det 1 /^' (D s t,m) and det (D ov ,m) 

We have seen, both in 2D and in 4D, that A^-fold nearly degenerate staggered eigenvalues 
manage to "mimic" a (single) overlap eigenvalue, e.g. ~ yjf with yjf = i(— ^at-i^fc) 1 ^ 2 m 2D. 
Obviously, such a relation cannot extend to the entire spectrum (they are rather different in the UV, 
and the number of modes does not match), but the relevant question is whether it would extend to 
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Figure 14: Scatter plot of j logdet(D st m ) versus logdet(D ov ,m) in a box of fixed physical volume and with 
m/g = 0.1 fixed. The first column is without filtering, the second with. For comparison the correlation of 
the overlap determinant with and without filtering is shown and ditto for the rooted staggered determinant, 
where the deviations are known to be 0{a 2 ) effects. The black lines represent the identity. Figure from 



determinant ratios (from a configuration U to U') as seen by these formulations, i.e. whether 

X\ %2 ••• 

A j A> ... 

If such a relation can be shown, specifically if one can prove that the deviation from an exact 
identity is caused by 0(a 2 ) cut-off effects, one would know that the rooted staggered determinant 
creates the "same" ensemble as the overlap action. 

In 2D it is rather easy to compute exact determinants, in particular if one has all eigenvalues on 
several ensembles in a matched physical volume on disk. This being the case with our Schwinger 
model studies, we included a correlation plot in |40| ], see Fig. [l4|. The sign-flipped contribution 
to the effective action on typical quenched backgrounds, as seen by (rooted) staggered or (single) 
overlap fermions is plotted against each other. In either case the determinant has been normalized 
by an artificial reference configuration which realizes the ensemble mean in both directions, and 
the physical box volume and quark mass are the same in all graphs. On coarse lattices, there 
is a rather loose correlation, i.e. we see a "cigar" shaped cloud. The interesting observation is 
that the thickness of the "cigar" shrinks on finer lattices, i.e. the correlation improves towards the 
continuum. This is particularly obvious if the two operators are filtered (second column). The 
full line always represents the identity, not a fit. To interpret this plot it is useful to recall that 
a correlation between So and — logdet(D ov . m ) or — ^logdet(.D st)m ) looks rather similar on coarse 



7i 72 



(7.1) 
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lattices, but here the correlation does not improve towards the continuum 2 . For comparison the 
third and fourth columns correlate different filtering levels for the overlap or staggered action, 
respectively. The point is that in these cases the difference is known to be an 0(a 2 ) effect. It is 
amusing to see that a filtered overlap determinant creates an ensemble that is closer to the one 
from a rooted staggered determinant at the same filtering level than to the one from an unfiltered 
overlap determinant. In summary, the similarity between the second and third column makes it 
plausible that det 1//iV ' (D st , m ) and det(D ov m ) create ensembles that differ by cut-off effects only (a 



more quantitative analysis is included in [40]). To my knowledge Fig. 14 contains the strongest 
evidence available right now that the overlap operator D ov might actually be (one version of) the 
desired "candidate" action D ca m , i.e. that the conjectured identity 

det(D 0V ,,„) = const • det 1 ^' (D stjm ) (l+<9(a 2 )) (7.2) 

might hold true with N t = 2 d l 2 (it is understood that the two operators have the same renormalized 
mass). On the other hand, it could be that all this is just misleading. In that case the only lesson to 
be learned would be that numerical evidence - no matter how convincing it looks - is never good 
enough to decide on conceptual issues. 

8. Staggered Chiral Perturbation Theory 

The last piece of evidence in favor of the rooting procedure comes from fitting hadron data to 
functional forms specific to the the low-energy effective theory for staggered quarks. For a general 



review on effective theories at finite lattice spacing I refer to [48] 



Let me first recall some of the foundations of standard Chiral Perturbation Theory (XPT) as 



formulated by Gasser and Leutwyler [49, 50]. The key observation is that in QCD at low energy 



(say below lGeV) chiral symmetry seriously constrains the interaction of "pions" (K,K,ri) with 



each other and with heavier degrees of freedom (e.g. nuclei, see the talk by U.-G. Meissner [51]). 



The way how this is implemented is through the Goldstone theorem [|52|, |53J]. More specifically, 
they start from the fundamental assumption that in QCD with Nf light quarks (the focus is on the 
meson sector with Nf = 2, 3) the pattern of spontaneous symmetry breaking (SSB) is 

SU(N f ) L x SU R (N f ) R — ► SU(N f ) v (8.1) 

and build the most general effective Lagrangian consistent with this (and a few more constraints 
like Lorentz covariance and parity conservation). In this framework an amplitude like, e.g. the 1=0 
Tin forward scattering amplitude [s = 0, v = (t—u)/(4M 7[ )] has a low-momentum expansion 

M 2 

Fnn{v) = T I=o {0,2M K {M K +v),2M K {M K -v)) = — -f + 0(p 4 ) (8.2) 



in powers of p . As announced above, ( |8.2P tends to zero for M%, V — > 0. The counting rule p ~m 
(where m is identified with M 2 ) is dictated by the reaction to infinitesimally small breaking terms. 

Given this background, it is clear that an extension of XPT to staggered quarks must build on 
an analysis how the broken taste symmetry at finite lattice spacing affects the pion mass. This was 



2 It simply means that there is a genuine unquenching effect that cannot be mimicked through a shift in j3. 
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done in a paper by Lee and Sharpe [54], and Fig. 15 illustrates their main findings. The generalized 
"pions" built, say, from an up-quark and a down-antiquark fall into different taste irreps d{ys®T:)u. 
Only the taste pseudoscalar (t = Ts, "P" in Fig. 15) is massless in the chiral limit at finite lattice 
spacing; this state is to be identified with the physical 7i + . All taste non-P states ("A,T,V,S" in 
Fig. 15) stay massive in the chiral limit; the excess of their squared masses over the z = T5 partner 



is nearly independent of the quark mass. Physical correlators project on taste pseudoscalars in 
the asymptotic states, but in intermediate loops all taste irreps do, of course, contribute. Even 
with intermediate flips into unphysical tastes being unavoidable, the continuum limit might be OK, 
if the contribution from physical intermediate states is far less than what continuum XPT would 



predict. The non-P mass excess diminishes like a 2 (see Fig. 15). Hence a scheme in which a cut 
in a 4-point function is built up from several contributions, each of which has a fraction of the 
continuum weight, would be OK if the fractions add up to 1. In fact, such a desired behavior is 
built into the effective theory via an enlarged symmetry group [and a dedicated prescription]. The 
key assumption is as follows. With Nf flavors of (4-taste) quark fields the pattern of SSB is 



SU{AN f ) L xSU{ANf) R -» SU(4Nf)\ 



(8.3) 



leading to 16N/ 2 — 1 pseudo-Goldstone bosons, collected in the 12 x 12 matrix (for Nj = 3) 



U 



where 



f<t> u n + K+\ 9jl6 
71- <D d K° = £ & h ^-T b 



(8.4) 



that transforms as U — > V L UV# under chiral rotations with unitary Vl,Vr. 

Staggered Chiral Perturbation Theory (SXPT) is the systematic extension of the Lee-Sharpe 



analysis to higher loop order, worked out in detail by Aubin and Bernard [55, 56]. The counting 
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rule p 2 ~ m ~ a 2 follows from the discussion above and the leading order Lagrangian reads 

-tv{MU+MU f ) + — 



f 2 £ 2m 2 

L = J—tr^Ud^) - -tr(MC/+Mf/ t ) + -^(<Z> u TS+ <Z> dJS+ <t> sJS ) + a 2 V TB (8.5) 



with /~ 122MeV the pion decay constant in the chiral limit, £~ (270 MeV) 3 the chiral condensate 
[in (MS,2GeV) conventions] and M = diag(m„, ...,/n^, ...,m s , ...) the 12 x 12 quark mass matrix. 



From a look at ( |8.4| , |8.5D it is clear that some details of the theory are rather different from ordinary 



XPT. The flavor sum in (^4) extends up to 9, not just 8 (here, X = ^ s /2/3^). This is necessary, 
since (jDl) yields 12 2 - 1 = 143 Goldstone and non-Goldstone bosons. With only 8T6 = 128 terms 
the dimensionality of the tangent space would be too small. On the other hand, with 9- 16 = 144 
terms there is space for one extra degree of freedom. This extra degree, the physical 77', gets lifted 
through the piece with the m 2 , prefactor in (3.5) which is designed to affect only the combined 



flavor and taste singlet state. That the T]' is an explicit degree of freedom creates a similarity with 
quenched or partially quenched XPT, even for ra se a = «?vai- The taste breaking potential Vtb in (fOb 



implements the excess mass of the non-Goldstone pions that has been illustrated in Fig. 15. With 



this setup several important quantities have been calculated at 1-loop level, e.g. for M n ,f K ,MK,fK 



it is known how they react to changes in m and a [55, 56]. 



The SSB pattern ( |§31 ) has been declared an assumption. Here, I would like to emphasize that 



the analogous pattern ( |8.1| ) remains an assumption in ordinary XPT, too. Of course, Gasser and 
Leutwyler mention that there is a wealth of evidence [from pions in the real world being so light, 
and even from the lattice] that the breaking ( |8.1| ) is actually realized in nature. The point is that their 
argument is evidence based, there is no proof that standard XPT is the correct low-energy theory 
of QCD. With SXPT the situation is rather analogous. One can compare the predictions to lattice 
data and see whether they fit. The MILC collaboration has done this in great detail, performing 



correlated fits in various circumstances (see e.g. [57]). Given the large number of free parameters, it 
is not so surprising that they get horrific covariance matrices, but in terms of the # 2 /d.o.f. obtained 
their fits are quite acceptable. Perhaps the most significant observation is that they cannot fit the 
taste pseudoscalar part of their data with continuum XPT, it is crucial to analyze their data in the 



full SXPT framework. Another point covered in [ 57 ] is that they generalize their formulae to allow 



for an arbitrary number of remaining tastes per flavor (not just 1). Fitting some of their lighter 
data with this extra parameter, they obtain 1.28(12) as the preferred number. More detailed tests 
have been proposed by Sharpe and van de Water and experiences in more exotic sectors have 



been reported at this conference [59, 60]. At this point there is a huge body of evidence that SXPT 
allows for an excellent description of a wealth of lattice data with the rooted staggered determinant. 
In principle there is, of course, the possibility that the "rooting trick" is incorrect and SXPT just a 
stunningly accurate model. Personally, I do not consider any accidental agreement of two flawed 
theories very likely, and this is why I rate the overall outcome of the MILC fits good evidence in 
favor of both the "rooting trick" and the SSB pattern (|8j|). Still, the generic comment applies that 
even the best numerical evidence does not amount to a mathematical proof. 

9. Summary 

The question whether in QCD with Nf = 2+1 light quarks the staggered action yields the right 
continuum limit remains unsolved. Instead of drawing any conclusion, let me just summarize: 
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1. Full QCD with =2+1 staggered fermions is controversial, since the Boltzmann weight 
det 1//4 (D st ) assumes a taste symmetry which is only approximate. 

2. Formally, the taste symmetry breaking is due to irrelevant interaction terms and should thus 
go away in the continuum limit. 

3. Weak coupling, filtering, RG blocking reduce the taste splitting and give staggered quarks 
more appealing features, but there is no guarantee that no trace is left in the continuum. 

4. One legal 1-flavor D ca m operator with det(D ca m ) = const •det 1//4 (D st m ) (\+0{a 2 )) is suffi- 
cient to re-interpret existing staggered ensembles as being generated with a local action. 

5. Even if such a D ca „, exists, the problem of (exact) unitarity in the fundamental theory remains 
when sticking to staggered spectroscopy. If the construction is exact, i.e. if there is no 0(a 2 ) 
term above, the problem is easily overcome by using this D ca m in the valence sector. 
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